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Abstract  This  paper  presents  the  development  of  a 
dynamic  calibration  technique  for  thermal  shear-stress 
sensors  using  acoustic  plane  wave  excitation.  The 
technique  permits  the  independent  variation  in  the 
mean  and  fluctuating  shear  stresses.  The  theoretical 
development  and  the  practical  implementation  of  the 
technique  are  presented.  The  studied  configuration  has 
the  capability  to  dynamically  calibrate  shear-stress 
sensors  up  to  20  kHz.  An  illustrative  application  of  this 
technique  to  an  uncompensated  silicon  micromachined 
thermal  shear-stress  sensor  operated  in  constant  current 
mode  is  discussed.  Specifically,  the  sensor  has  been 
statically  calibrated  over  a  range  of  wall  shear  stress 
from  7  to  80  mPa.  A  dynamic  calibration  of  the  sensor 
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over  a  range  of  2-12  mPa  has  been  performed  up  to 
7  kHz. 


1  Introduction 


The  measurement  of  wall  shear  stress  in  a  turbulent  flow 
provides  insight  into  complex  flow  phenomena,  including 
viscous  drag,  flow  separation/reattachment  and  transi¬ 
tion  to  turbulence.  At  high  Reynolds  numbers,  the  spatial 
length  scale  of  the  fluid  structures  of  interest  can  be  on  the 
order  of  hundreds  of  microns  or  less,  necessitating  a  sen¬ 
sor  element  size  of  the  same  order  to  minimize  spatial 
averaging.  Also,  to  accurately  capture  the  complete 
spectrum  of  turbulent  fluctuations,  a  shear-stress  sen¬ 
sor  with  usable  bandwidth  of  over  1  kHz  is  required. 
For  example,  it  is  well  known  that  the  ratio  of  the 
boundary  layer  thickness  3  to  the  Kolmogorov  microscale 
rj  scales  as 


n 

3 


Re 


-3/4 


(i) 


where  Res  =  u3/v  is  the  Reynolds  number,  u  is  the  typical 
eddy  velocity  scale  (u/U~  10“  2),  U  is  the  boundary  layer 
edge  velocity,  and  v  is  the  kinematic  viscosity  (Tennekes 
and  Lumley  1972).  Similarly,  the  ratio  of  the  convective 
time  scale  to  the  Kolmogorov  time  scale  T  scales  as 


(2) 


For  a  relatively  low  Reynolds  number,  low-speed  tur¬ 
bulent  boundary  layer  possessing  an  edge  velocity  of 
(7=50  m/s  and  a  thickness  3  =  5  cm,  the  Kolmogorov 
scales  are  rj=  100  |im  and  T=  2  ms.  The  stringent  spatial 
and  temporal  resolution  requirements  naturally  point 
towards  microfabricated  sensors  as  a  means  of  achieving 
the  required  level  of  performance. 

Conventional  as  well  as  micromachined  shear  stress 
measurement  techniques  can  be  broadly  classified  into 
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direct  and  indirect  methods.  Several  microfabricated 
shear-stress  sensors  of  both  the  direct  (Padmanabhan 
et  al.  1996,  1997;  Pan  et  al.  1999)  and  indirect  (Oudh- 
eusden  and  Huijsing  1988;  Liu  et  al.  1994,  1999;  Cain 
et  al.  2000)  variety  have  been  reported.  The  direct  shear- 
stress  sensors  are  tethered  floating-element  devices  that 
employ  piezoresistive  (Ng  1990;  Goldberg  et  al.  1994), 
optical  (Padmanabhan  et  al.  1996,  1997)  and  capacitive 
(Schmidt  et  al.  1988;  Pan  et  al.  1999)  transduction 
mechanisms.  Indirect  shear-stress  sensors  are  tempera¬ 
ture-resistance  transducers  that  operate  on  heat  transfer 
analogies.  Recently,  Lofdahl  and  Gad-el-Hak  (1999) 
presented  a  review  of  MEMS  pressure  and  shear-stress 
sensors  for  turbulent  flows.  Naughton  and  Sheplak 
(2002)  expanded  that  contribution  by  providing  an  up¬ 
dated  list  of  references  and  discussing  some  of  the  cur¬ 
rent  misconceptions  and  controversies  involving  the  use 
of  MEMS-based  shear-stress  sensors  for  quantitative 
measurements. 

To  accurately  capture  the  complete  spectrum  of  shear 
stress  in  a  turbulent  flow,  a  sensor  should  possess  a 
known  frequency  response  function.  An  ideal  shear 
stress  measurement  system  would  possess  equal  static 
and  dynamic  sensitivities  with  a  flat,  minimum-phase 
frequency  response  over  all  frequencies.  In  practice, 
however,  the  inherent  compliance,  inertance  and  dissi¬ 
pation  of  the  system  limit  the  useful  bandwidth  of  the 
sensor.  For  example,  the  usable  bandwidth  of  a  thermal 
shear-stress  sensor  is  determined  by  the  thermal  inertia 
of  the  substrate  (Bellhouse  and  Schultz  1967,  1968). 
Bellhouse  and  Schultz  describe  how  the  dynamic  re¬ 
sponse  is  further  complicated  by  the  frequency  depen¬ 
dent  heat  conduction  into  the  substrate.  This  creates  a 
low-frequency  roll-off  in  the  gain  factor  of  the  frequency 
response  function,  as  well  as  a  frequency  dependent 
phase  lag,  resulting  in  unequal  static  and  dynamic  sen¬ 
sitivities.  The  low-frequency  substrate  conduction  losses 
also  limit  electronic  signal  injection,  commonly  used  in 
hot-wire  anemometry,  to  providing  qualitative  infor¬ 
mation.  Since  the  exact  character  of  the  system  is  un¬ 
known,  the  response  of  such  a  system  to  electronic  signal 
injection  is  unclear.  Consequently,  electronic  signal 
injection  fails  to  yield  an  accurate  estimate  of  the  fre¬ 
quency  response  function  (Freymuth  1981).  Therefore, 
to  minimize  uncertainties  in  correlation  and  spectral 
analysis,  an  in  situ  dynamic  calibration  technique  is  re¬ 
quired  to  determine  the  frequency  response  function. 

A  dynamic  calibration  technique  for  thermal  shear- 
stress  sensors  requires  a  known  mean  and  broadband 
shear  stress  input  which,  unfortunately,  is  difficult  to 
realize  in  practice.  Bellhouse  and  Schultz  (1968)  devel¬ 
oped  a  dynamic  calibration  technique  based  on  the 
streamwise  oscillation  of  a  small,  flat  plate  containing  a 
flush  mounted  sensor,  thus  effectively  modulating  the 
freestream  velocity.  This  technique  permitted  a  dynamic 
calibration  range  of  200-1200  Hz.  Chew  et  al.  (1998) 
have  designed  an  experimental  set  up  that  uses  Couette 
flow  between  a  recessed  rotating  disc  and  a  flat  sta¬ 
tionary  disc  to  characterize  the  dynamic  response  of  wall 


mounted  hot-wire  and  hot-film  probes.  This  technique 
was  used  to  dynamically  calibrate  flush  mounted 
hot-film  probes  up  to  2.5  kHz.  One  limitation  to  this 
calibration  technique  is  that  the  mean  shear  stress  and 
frequency  of  the  disturbance  are  functions  of  the  rota¬ 
tional  speed  for  a  fixed  plate  geometry  and  thus  are  not 
independently  adjustable.  More  recently,  a  dynamic 
calibration  technique  utilizing  acoustic  waves  to  gener¬ 
ate  a  known  sinusoidal  wall  shear  stress  under  zero- 
mean  flow  conditions  has  been  reported  (Sheplak  et  al. 
2001).  This  technique  was  used  for  the  dynamic  cali¬ 
bration  of  a  micromachined  floating-element  sensor  that 
exhibits  a  linear  response  to  shear  stress.  This  procedure 
is  not  suitable  for  a  thermal  sensor  that  possesses  a  non¬ 
linear  response  and  requires  a  mean  flow  to  ensure 
forced  convection  of  heat  from  the  sensor.  This  paper 
presents  the  development  of  an  acoustic  calibration 
technique  with  independent  variation  in  the  mean  and 
fluctuating  shear  stresses  to  enable  static  and  dynamic 
calibration  of  thermal  shear-stress  sensors.  The  theo¬ 
retical  development  and  the  practical  implementation  of 
the  dynamic  calibration  technique  are  presented.  The 
outline  of  the  paper  is  as  follows.  The  next  section  pre¬ 
sents  the  theoretical  background  of  the  technique.  This 
is  followed  by  a  proof-of-concept  experiment  using  an 
uncompensated  silicon  micromachined  thermal  shear- 
stress  sensor  operated  in  the  constant  current  mode. 
Finally,  conclusions  are  drawn  and  suggestions  for 
future  work  are  given. 


2  Background 

The  response  of  the  hot  film  sensor  to  dynamic  shear 
fluctuations  is  band-limited  due  to  the  inherent  dissipa¬ 
tive  and  storage  mechanisms  of  the  system.  The  energy 
balance  of  a  hot  film  is  similar  to  that  of  a  conventional 
hot  wire  except  for  the  frequency-dependent  heat  con¬ 
duction  into  the  substrate  governed  by  the  temperature 
gradient  at  the  sensor/substrate  interface.  For  the  case  of 
hot  wires,  the  conduction  of  heat  into  the  supporting 
prongs  is  negligible  compared  to  the  convective  losses 
for  length-to-diameter  ratios  of  the  order  1,000  (Perry 
et  al.  1979).  However,  in  the  case  of  hot  films,  this 
conduction  term  is  always  significant  and  cannot  be 
accurately  modeled  analytically  due  to  the  complex 
nature  of  the  conjugated  heat-transfer  problem  (Comte- 
Bellot  1976).  This  is  even  more  complex  for  sensor  de¬ 
signs  possessing  the  thermal  sensing  element  on  a  thin 
diaphragm  that  is  stretched  over  a  vacuum  cavity  (Liu 
et  al.  1994,  1999;  Cain  et  al.  2000). 

It  is  known  that  thermal  shear-stress  sensors  exhibit  a 
non-linear  response  to  shear  stress.  If  the  system  can  be 
modeled  approximately  as  locally  linear  and  time 
invariant,  then  a  direct  dynamic  calibration  is  possible 
provided  that  a  known  fluctuating  shear  stress  input  can 
be  produced.  For  the  linearity  assumption  to  hold,  the 
shear  stress  perturbation  (T'waii)  must  be  sufficiently 
smaller  than  the  established  mean  shear  stress  (f),  such 
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that  V[  =  (dF/df)T^a//,  where  V'T  is  the  corresponding 
voltage  fluctuation  due  to  the  fluctuating  shear  stress  and 
dF/df  is  the  local  static  sensitivity.  In  the  following 
sections,  the  theoretical  basis  for  the  dynamic  calibration 
using  acoustic  excitation  in  a  square  duct  is  presented. 


2.1  Stokes  layer  excitation 


The  basic  principle  of  this  technique  relies  on  the  fact 
that  the  particle  velocity  of  the  acoustic  waves  must  be 
zero  at  the  wall  due  to  the  no-slip  boundary  condition. 
This  leads  to  the  generation  of  a  frequency-dependent 
boundary  layer  thickness  (5(<x>)  and  a  corresponding  wall 
shear  stress  T'wan(co).  The  non-dimensional  thickness  of 
this  viscous  region,  known  as  the  Stokes  layer,  is  given 
by  Panton  (1996)  as 


where  co  is  the  angular  frequency,  b  is  an  appropriate 
length  scale,  and  r\  =  y/cob2/v  is  the  non-dimensional 
Stokes  number. 

Now  consider  a  set  of  plane,  traveling  acoustic  waves 
bounded  by  the  walls  of  a  rectangular  duct  of  cross- 
section  2ax2b  with  (a>b),  as  shown  in  Fig.  1.  The  waves 
travel  in  the  positive  v-direction  and  induce  a  velocity 
perturbation  in  the  fluid.  For  the  case  of  purely  traveling 
acoustic  plane  waves  superimposed  on  a  fully  developed 
mean  flow  with  centerline  axial  velocity  U<^  c,  where  c  is 
the  isentropic  speed  of  sound,  the  linearized  compress¬ 
ible  momentum  equation  in  the  axial  direction  reduces 
to  the  classic  problem  of  a  duct  flow  driven  by  an 
oscillatory  pressure  gradient  (Sheplak  et  al.  2001).  The 
unsteady  momentum  equation  in  differential  form  is 
given  by 


du 

~di 


1  dp  2 
-^-  + vV2m, 

p  OX 


(4) 


where  p  is  the  density  of  the  fluid,  u  the  axial  velocity, 
and  p  is  the  pressure.  We  assume  that  the  acoustic  wave 
is  a  traveling  plane  wave  in  the  +  v  direction 


p{x ,  t)  =  p  Re 


^((ot—kx—n/2) 


(5) 


where  p'  is  the  amplitude  of  the  pressure  perturbation, 
Re  [  ]  denotes  the  “real part  of[  ]”,  j  =  and  k  =  colc 
is  the  acoustic  wavenumber.  The  assumed  form  in  Eq.  5 
translates  to  an  oscillatory  pressure  gradient  of  the  form 


— =  K cos  (cot  —  kx ), 
p  dx 


(6) 


where  K=p'co/pc.  The  velocity  fluctuation  is  then  of  the 
form 


(x,y,z,t)  = 


(7) 


where  u  (y,  z)  is  the  complex  amplitude  of  the  velocity 
perturbation.  The  velocity  fluctuation  is  subject  to  no¬ 
slip  boundary  conditions 

u  =  0  at  z  =  ±a  and  u  =  0  at  y  =  ±b.  (8) 


Solutions  to  this  problem  have  been  obtained  by 
Drake  (1965),  Fan  and  Chao  (1965),  and  O’Brien  (1975). 
For  example,  the  fluctuating  velocity  on  the  central 
plane  of  the  duct  (z  =  0)  can  be  obtained  from  O’Brien’s 
solution  as  the  real  part  of 


'(v,o,0=— ei(“'“fa“(7t/2))x 

pc 


(-i  r 

Pn 


COS (Pny/b)  |  cosh (j^y/b) 
cosh  (yna/b)  cosh  (pn) 


(9) 


where  p„  =  ((2n+l)/2)n  ,  qn  =  ((2n  +  1  )j2)n{bja),  yn  = 

Vrf+Pf,  and p„  =  ^/ql+jri2.  As  a/b  ->  the  solu- 
tion  tends  to  the  well-known  result  of  flow  in  a  one¬ 
dimensional  slot  subject  to  an  oscillatory  pressure-gra¬ 
dient 


0  =  1 1  _ cosh  M. b-Ml }.  (10) 

pc  {  cosh  (riy/j)  J 

Figure  2  shows  a  plot  of  the  normalized  disturbance 
velocity  magnitude  along  the  centerline  versus  y/b  for 
various  Stokes  numbers  (rj)  in  a  square  ( a  =  b)  duct.  Both 
the  one-dimensional  and  two-dimensional  solutions  are 
shown.  Note  that,  for  square  ducts,  there  is  a  significant 
difference  between  the  one-dimensional  and  two- 
dimensional  solutions  only  for  small  Stokes  numbers, 
which  correspond  to  low-frequencies.  As  the  Stokes 
number  increases,  the  decreasing  thickness  of  the  Stokes 
layer  becomes  apparent,  in  accordance  with  Eq.  3. 

It  is  important  to  note  that  Fan  and  Chao  provide 
only  the  real  part  of  the  solution,  which  is  insufficient  to 
calculate  the  complex  shear  stress.  In  addition,  by 
comparing  the  various  solutions,  the  series  solution  in 
Drake  was  found  to  be  inaccurate  at  high  Stokes  num¬ 
bers,  while  the  series  solution  provided  by  O’Brien  does 
not  converge  at  low  Stokes  numbers  ( r\  <  10). 


Fig.  1  Illustration  of  the  acoustic  waveguide  and  coordinate  axes 
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Fig.  2  Normalized  magnitude  of  the  velocity  profile  for  various 
Stokes  numbers  (rj)  for  one-dimensional  and  two-dimensional 
oscillatory  rectangular  duct  flow 


The  fluctuating  wall  shear  stress  induced  by  the 
oscillatory  pressure  gradient  at  the  center  of  the  lower 
wall  of  the  duct  is  given  by 


^walfiD 


dd_ 

dy 


y=—b,z=0 


(ii) 


Using  the  two-dimensional  solution  of  O’Brien,  Eq.  9, 
the  complex  fluctuating  shear  stress  is 


Twall2D 


—2vp' 

be 


oo 

qH cot-kx-n/2 ) 

n— 0 


(-1)” 

Pn 


Pn  sin  (p„) 
cosh  ( yna/b ) 


Pn  tanh  (Pn) 


(12) 


Since  the  one-dimensional  and  two-dimensional  solu¬ 
tions  are  essentially  identical  for  rj  >  2  (which  as  shown 


later  corresponds  to  the  test  conditions  in  the  current 
study),  the  one-dimensional  fluctuating  wall  shear  stress 
is  used 


^  wall  id 


Cm) 


(13) 


The  above  equation  indicates  that  the  amplitude  of  the 
shear  stress  is  directly  proportional  to  the  product  of  the 
pressure  perturbation  and  the  square  root  of  the  exci¬ 
tation  frequency.  A  simple  algebraic  manipulation 
shows  that  the  fluctuating  shear  stress  amplitude  is 
proportional  to  the  Stokes  number  since 


OJV 


[ip 

b  pc 


(14) 


Figure  3  shows  the  magnitude  |T'wan///|  and  phase  of  the 
fluctuating  shear  stress  with  respect  to  the  pressure 
perturbation  for  the  one-dimensional  case.  The  linear 
dependence  versus  Stokes  number  for  rj>  1  is  revealed  in 
the  figure.  Also  note  that  for  rj  >  4,  the  shear  stress  leads 
the  pressure  by  45°.  Equation  13  is  valid  for  plane 
progressive  waves.  At  frequencies  higher  than  the  cut-on 
frequency  of  the  acoustic  waveguide,  non-planar  modes 
propagate  within  the  duct,  in  addition  to  plane  waves 
(Morse  and  Ingard  1968).  Therefore,  the  bandwidth 
over  which  the  solution  is  valid  is  limited  by  the  wave¬ 
guide  characteristics  of  the  tube.  Furthermore,  the 
acoustic  waveguide  must  be  properly  terminated  to 
eliminate  reflected  waves.  Details  of  the  experimental 
setup  are  provided  below. 


2.2  Plane  wave  propagation  in  a  square  duct 

Wave  propagation  in  a  duct  is  governed  by  the  solution 
to  the  wave  equation.  For  one-dimensional  wave  prop¬ 
agation,  the  width  of  the  duct  must  be  small  in 


Fig.  3  Magnitude  and  phase  of 
the  ratio  of  the  fluctuating  wall 
shear  stress  with  respect  to  the 
pressure  perturbation  versus 
Stokes  number  (rj)  for  one¬ 
dimensional  oscillatory 
rectangular  duct  flow 
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comparison  to  the  acoustic  wavelength.  This  can  be 
verified  by  solving  the  wave  equation  in  the  duct 
(Dowling  and  Williams  1983) 


1  d2p' 
c 2  dt 2 


-  vV  =  0. 


(15) 


Consider  the  case  of  a  square  duct  with  rigid  walls 
and  cross-sectional  length  2a.  The  solution  to  the  wave 
equation  for  an  acoustic  wave  of  frequency  co  for  the 
( m,n )  mode  is 


p\x,y,z,t)  =  cos 


mn(y  +  a) 
2  a 


cos 


nn(z  +  a) 
2  a 


A  J(cot-kmnx )  D  J((Dt+kmnx) 


(16) 


The  corresponding  wall  shear  stress  (y  =  a)  is 


_  8a dp  ^  /_iyj 

1)_ 


X 


cosh  ( nnz/2a ) 
cosh(?77r/2) 


(—  sin  {nn/2)) 


(20) 


Figure  4  shows  the  variation  of  shear  stress  along  one  of 
the  walls  of  the  square  duct.  It  is  seen  that  the  maximum 
shear  stress  occurs  at  the  center  of  the  wall  with  only  a 
1%  decrease  for  1  mm  deviation  from  the  centerline  for 
an  8.5  mm  square  duct.  Compared  to  the  sensing  ele¬ 
ment  length  of  0.2  mm,  we  conclude  that  the  spanwise 
shear  stress  variation  across  the  sensor  is  negligible. 


where  kmn  is  the  wave  number  given  by  the  dispersion 
relation  (Dowling  and  Williams  1983) 


When  kmn  is  real,  Eq.  16  represents  traveling  waves 
propagating  in  the  x-direction.  If  kmn  is  imaginary,  it 
represents  non-propagating,  evanescent  waves  whose 
magnitudes  decay  exponentially  with  axial  distance. 
The  fundamental  mode  (0,0)  is  a  plane  wave  that  will 
always  propagate,  but  the  higher  order  modes  will  be¬ 
come  non-evanescent  only  when  the  cross-sectional 
width  of  the  duct  is  greater  than  half  of  an  acoustic 
wavelength,  i.e.,  2a  >2/2.  For  example,  the  cut-on  fre¬ 
quency  for  the  first  higher  order  modes  of  an  8.5  mm 
square  tube  is  20  kHz. 

As  mentioned  earlier,  the  Stokes  layer  excitation  is 
superimposed  on  an  established  mean  flow.  The  static 
sensitivity  of  the  sensor  limits  the  magnitude  of  the  shear 
stress  perturbation  that  can  be  applied  for  the  local 
linear  assumption  to  hold.  The  determination  of  mean 
wall  shear  stress  and  the  optimal  sensor  placement 
within  a  square  duct  are  discussed  in  the  remainder  of 
this  section. 


3  Experimental  set-up 

The  experiments  were  carried  out  in  the  Interdisciplinary 
Microsystems  Laboratory  at  the  University  of  Florida. 
A  plane  wave  tube  (PWT)  capable  of  generating  purely 
traveling,  acoustic  plane  waves  up  to  20  kHz  was  de¬ 
signed  and  fabricated.  Two  series  of  experiments  were 
performed,  one  to  characterize  the  response  of  the  PWT, 
the  other  to  demonstrate  the  Stokes  layer  dynamic  cal¬ 
ibration  technique  for  thermal  shear-stress  sensors. 

Dynamic  calibration  of  the  shear-stress  sensor  first 
requires  a  thorough  static  characterization  that  includes 
current-voltage,  resistance-temperature  and  static  wall 
shear  stress  sensitivity  experiments.  A  detailed  descrip¬ 
tion  of  these  experiments  and  related  equipment  is  given 
by  Cain  (1999),  Cain  et  al.  (2000),  and  Sheplak  et  al. 
(2002). 


3.1  Plane  wave  tube 

A  schematic  of  the  acoustic  waveguide  setup  is  shown 
in  Fig.  5.  The  aluminum  tube  has  a  length  of  1.45  m,  a 
square  duct  cross-section  of  8.5  mm,  and  a  wall  thick¬ 
ness  of  1.27  cm.  Two  JBL  2426J-compression  drivers 


2.3  Laminar  flow  in  a  square  duct 


For  an  incompressible,  steady,  laminar,  fully  developed 
flow  in  the  square  duct,  the  axial  momentum  conserva¬ 
tion  equation  reduces  to  Poisson’s  equation 


1  dp  d2u  d2u 
fi  dx  dy2  +  dz 2 


The  solution  to  the  above  equation  for  the  velocity 
distribution  in  the  duct  is  (White  1991) 
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Fig.  4  Mean  shear  stress  variation  along  a  wall  of  the  square  duct 
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Fig.  5  Schematic  of  the  plane  wave  tube  setup 


are  attached  to  the  plane  wave  section  via  transition 
ducts  with  varying  inner  cross-section.  The  1/8”  Bruel& 
Kjaer  (B&K,  Denmark)  4138  microphone  and  thermal 
shear-stress  sensor  are  flush  mounted  on  opposite  sides 
of  the  tube  wall  at  the  same  axial  distance  (1.0  m)  from 
the  speakers.  The  B&K  4138  is  connected  to  a  B&K 
2669  preamplifier  and  powered  by  a  B&K  2804  power 
supply.  A  Stanford  Research  Systems  (Sunnyvale,  CA) 
SR785  spectrum  analyzer  serves  as  the  signal  source 
and  the  data  collection  unit.  The  speakers  are  driven  by 
a  Techron  (USA)  7540  two-channel  amplifier.  The  mi¬ 
crophone  and  the  sensor  are  housed  in  custom-ma¬ 
chined  aluminum  and  Lucite  packages,  respectively 
(Chandrasekaran  2000).  The  Lucite  package  (Fig.  6) 
consists  of  a  multi-step  circular  terrace  structure  with 
an  O-ring  to  provide  an  airtight  seal  and  copper  rods  to 
provide  front-to-back  electrical  contacts.  The  sensor  is 
bonded  to  the  package  using  fast-cure  epoxy  within  a 
precision-cut  cavity  such  that  flush  mounting  is 
achieved.  Electrical  contact  is  made  between  the  sensor 
and  copper  rods  using  an  ultrasonic  ball  wire-bonder 
with  0.025  mm  diameter  gold  wire.  When  the  packaged 
sensor  is  mounted  in  the  PWT,  the  wire  bonds  are 
oriented  downstream  of  the  active  area  of  the  sensor  to 
minimize  any  disruption  to  the  flow.  Both  ends  of  the 
PWT  terminate  into  compliant  tubing  to  approximate  a 


non-reflective  acoustic  termination  (Morse  and  Ingard 
1968).  The  tubing  is  a  duct  with  a  compliant  wall  and  is 
intended  to  dissipate  the  acoustic  energy  propagating 
through  it.  However,  any  slight  change  in  acoustic 
impedance  at  the  interface  of  the  rigid  and  compliant 
ducts  causes  scattering  or  reflection  of  the  incident 
acoustic  waves  (Morse  and  Ingard  1968).  In  addition, 
any  change  in  cross-section  of  the  compliant  tubing 
also  modifies  the  acoustic  impedance.  The  reflected 
waves  propagate  in  the  negative  x-direction  resulting  in 
a  standing  wave  field  in  the  PWT.  If  the  effect  due  to 
the  impedance  mismatch  is  strong  it  could  lead  to  er¬ 
rors  in  the  shear  stress  measurement,  which  assumes  a 
purely  traveling  acoustic  wave.  However,  by  measuring 
the  reflection  coefficient  at  the  interface  of  the  rigid/ 
compliant  wall  the  assumed  form  of  the  pressure  per¬ 
turbation  can  be  modified  to  include  the  reflected  wave 
component. 

The  termination  tubing  on  the  inlet  (speaker)  end  of 
the  PWT  also  serves  as  the  entrance  for  the  mean  flow. 
The  flow  rate  is  manually  adjusted  with  a  Tescom  (USA) 
0-103  kPa  pressure  regulator  that  is  supplied  by  a 
compressed  air  source.  The  transducers  are  located  at  a 
distance  of  120  channel  heights  from  the  inlet,  ensuring 
that  a  fully  developed  laminar  flow  is  established  in  the 
duct  (Fox  and  McDonald  1998).  A  computational 
study  was  conducted  using  commercial  CFD  software 
(FLUENT)  to  validate  this  assumption  even  for  the 
“worst-case  scenario”  in  which  large  secondary  flow 
disturbances  were  introduced  in  the  entrance  section. 
The  results  of  the  simulations  indicate  that  for  a  laminar 
flow  the  fully  developed  assumption  is  valid  for  suffi¬ 
ciently  large  distances  down  the  duct,  regardless  of  the 
inlet  condition.  While  this  distance  may  vary  from  case 
to  case,  for  the  experiments  performed,  the  sensor  was 
well  within  a  region  where  the  fully  developed  assump¬ 
tion  provided  an  adequate  model  for  wall  shear  stress 
and  matched  with  CFD  results  to  within  0.3%.  The 
termination  tubing  at  the  other  end  of  the  PWT  is  open 
to  the  atmosphere.  Pressure  taps  located  at  regular 
spacing  along  the  length  of  the  test  section  are  used  to 
ascertain  the  streamwise  pressure  distribution  in  the 
PWT  to  verify  the  fully  developed  incompressible  mean 
flow  assumption  (Chandrasekaran  2000).  The  axial 
pressure-gradient  is  measured  via  a  0-124.5  Pa  Heise 
(USA)  differential  pressure  transducer  and  is  used  to 
calculate  the  mean  shear  stress  induced  via  Eq.  20.  A  30- 
gauge  type  K  thermocouple  is  mounted  2.54  cm  down¬ 
stream  of  the  sensor  between  the  top  and  bottom 
channel  surfaces  to  measure  the  fluid  temperature.  The 
fluid  viscosity  and  density  were  then  calculated  using 
Sutherland’s  formula  (White  1991)  and  the  ideal  gas  law, 
respectively. 

The  limiting  condition  for  laminar  flow  in  the  PWT 
was  theoretically  and  experimentally  determined.  The 
duct  flow  is  theoretically  laminar  for  Re  <2300,  corre¬ 
sponding  to  a  wall  shear  stress  of  90  mPa.  The  transition 
from  laminar  to  turbulent  flow  was  investigated  by 
monitoring  the  output  voltage  of  the  shear-stress  sensor 


Fig.  6  Shear-stress  sensor  package 


62 


while  increasing  the  mean  flow  rate.  The  sensor  voltage 
began  to  exhibit  significant  random  fluctuations  for 
mean  wall  shear  stresses  >100  mPa.  Once  the  flow  be¬ 
comes  turbulent,  Eq.  19  is  no  longer  an  accurate  repre¬ 
sentation  of  the  mean  velocity  profile  in  the  duct. 


3.2  Thermal  shear-stress  sensor 


The  test  device  structure  consists  of  a  1500A 
-thickx4  |im-widex200  pm-long  platinum  sensing 
element  on  top  of  a  1500A  -thick  silicon-rich  silicon 
nitride  membrane.  The  sensing  element  is  isolated  from 
the  bulk  silicon  substrate  via  a  200  pm-diameter, 
10  pm-deep  vacuum  cavity  (Cain  et  al.  2000).  This 
structure  reduces  the  sensitivity-limiting  thermal  con¬ 
duction  into  the  silicon  substrate.  The  use  of  a  vacuum 
cavity  to  reduce  substrate  conduction  effects  was  first 
reported  by  Liu  et  al.  (1994).  Each  sensor  possesses  two 
gold  leads  at  each  end  of  the  sensing  element  to  permit 
4-point  probe  characterization  exclusive  of  the  effects  of 
the  biasing  circuitry.  The  details  of  the  device  design, 
fabrication,  and  characterization  are  given  by  Cain 
(1999),  Cain  et  al.  (2000),  and  Sheplak  et  al.  (2002).  A 
plan-view  scanning  electron  micrograph  (SEM)  of  the 
device  is  shown  in  Fig.  7,  and  a  cross-sectional  sche¬ 
matic  is  shown  in  Fig.  8.  The  sensing  element  is  resis- 
tively  heated  to  a  temperature  greater  than  the  gas 
temperature  defined  by  the  non-dimensional  thermal 
overheat  ratio 


OH  = 


(21) 


where  Ts  and  Tg  are  the  absolute  sensor  and  gas  tem¬ 
peratures,  respectively.  The  convection  of  heat  from  the 
sensor  is  related  to  the  wall  shear  stress  by  the  Reynolds 
analogy  (White  1991),  and  is  measured  by  monitoring 
changes  in  the  temperature  dependent  resistance  of  the 
sensing  element.  Constant  current  sensor  excitation  and 
voltage  measurements  are  performed  with  a  Keithley 
(USA)  2400  source-meter. 
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Fig.  7  Plan-view  SEM  of  the  active  area  of  the  shear-stress  sensor 
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Fig.  8  Cross-sectional  schematic  of  the  shear-stress  sensor 


4  Experimental  results 

Details  of  the  static  calibration  of  this  device  are  given  in 
Cain  (1999)  and  Sheplak  et  al.  (2002).  This  section 
summarizes  the  results  of  static  and  dynamic  wall  shear 
stress  calibrations  of  the  micromachined  thermal  shear- 
stress  sensor. 


4.1  Static  calibration 

The  non-linear  static  response  of  the  sensor  for  varying 
mean  shear  stress  (10-80  mPa)  is  shown  in  Fig.  9.  The 
magnitude  of  the  output  response  increases  with  both 
increasing  mean  shear  stress  and  overheat  ratio.  In 
addition,  the  degree  of  non-linearity  in  the  static  re¬ 
sponse  is  significantly  less  than  that  of  a  hotwire 
(Comte-Bellot  1976).  The  local  slope  of  the  curve  rep¬ 
resents  the  static  sensitivity  of  the  device  at  a  particular 
operating  and  mean  flow  condition.  As  shown  in 
Fig.  10,  the  static  sensitivity  increases  with  increasing 
overheat  ratio,  but  decreases  with  increasing  mean  shear 
stress.  In  the  following  section,  results  of  the  dynamic 
calibration  of  the  thermal  shear-stress  sensor  in  a  con¬ 
stant  current  mode  of  operation  are  summarized. 


4.2  Dynamic  calibration 

In  our  experiment,  the  fluctuating  input  shear  stress  is 
determined  from  the  measured  acoustic  pressure  and 
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Fig.  9  Static  response  of  the  thermal  shear-stress  sensor  for  varying 
overheat  ratios 
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Fig.  10  Static  sensitivity  of  the  thermal  shear-stress  sensor  as  a 
function  of  mean  shear  stress 
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frequency  as  described  by  Eq.  13.  The  measured  fre¬ 
quency  response  function  is  normalized  by  the  static 
sensitivity,  dV jdx  to  yield  the  normalized  frequency  re¬ 
sponse  function  (Sheplak  et  al.  2001) 


H(co)  = 


V(co)  df 


(22) 


Since  the  magnitude  of  the  shear  stress  is  proportional  to 
the  product  of  the  acoustic  pressure  and  the  square  root 
of  the  excitation  frequency,  a  calibration  of  the  com¬ 
pression  driver  within  the  PWT  was  performed.  Fig¬ 
ure  11  shows  a  plot  of  the  magnitude  of  the  frequency 
response  of  the  speaker/PWT  system  for  a  constant 
voltage  excitation.  To  impose  a  known  shear  stress  on 
the  sensor  and  avoid  acoustic  non-linearities,  the  mag¬ 
nitude  of  the  speaker  driving  voltage  is  adjusted  to  ob¬ 
tain  a  nearly  constant  input  pressure  perturbation.  The 
oscillations  in  sound  pressure  level  towards  the  begin¬ 
ning  of  the  frequency  sweep  are  due  to  reflections  from 
the  ends  of  the  PWT.  Although  provisions  were  made  to 
eliminate  reflections,  the  length  of  compliant  tubing 
(3.5  m)  was  insufficient,  especially  at  lower  frequencies 
(<200  Hz). 


Frequency  (Hz) 


The  sensor  was  dynamically  calibrated  using  a  con¬ 
stant  amplitude  tone  (105  dB  ref  20  pPa),  at  an  overheat 
ratio  of  0.92  and  a  mean  shear  stress  of  50  mPa.  Fig¬ 
ure  12  shows  the  gain  factor  of  the  sensor  as  a  function 
of  frequency.  The  uncompensated  sensor  displays 
appreciable  sensitivity  up  to  4  kHz,  with  the  —  3  dB 
point  at  ~600  Hz.  Beyond  7  kHz,  the  response  is  on  the 
order  of  the  sensor  noise  floor,  which  in  this  case  indi¬ 
cates  the  baseline  voltage  output  of  the  sensor  before  the 
introduction  of  the  acoustic  waves.  A  sound  pressure 
level  of  105  dB  provides  theoretical  fluctuating  shear 
stress  amplitudes  of  approximately  2  mPa  at  200  Hz  to 
12  mPa  at  7  kHz.  For  these  levels,  the  local  linearization 
assumption  is  satisfied  (see  Fig.  9).  To  verify  the  linear 
relationship  between  the  output  voltage  fluctuation  and 
the  shear  stress  perturbation,  the  sensor  was  excited 
using  acoustic  waves  of  varying  amplitude  (90-120  dB) 
at  constant  frequencies.  Figure  13  confirms  this  linear 
relationship. 

Returning  back  to  Fig.  12,  the  static  and  dynamic 
response  should  ideally  be  equal,  resulting  in  a  constant 
gain  factor  of  OdB.  As  previously  mentioned,  the  ther¬ 
mal  inertia  of  the  shear-stress  sensor  will  cause  a  roll  off 
in  the  frequency  response  function.  In  Fig.  12,  the  roll¬ 
off  in  the  gain  factor  appears  to  be  ~40  dB/decade, 
indicative  of  a  highly  damped  second-order  system.  Fing 
and  Hubbard  (1956)  used  a  1 /2-order  system  to  model 
the  response  of  a  sensor  on  a  semi-infinite  medium.  The 
difference  in  the  response  of  micromachined  sensor  and 
the  conventional  sensor  described  in  Fing  and  Hubbard 
(1956)  can  be  explained  by  the  presence  of  a  sealed 
vacuum  cavity.  The  sealed  cavity  drastically  reduces  the 
unsteady  heat  conduction  into  the  substrate  because  the 
conduction  losses  are  limited  to  the  thin  nitride  mem¬ 
brane  as  opposed  to  the  entire  substrate. 

The  phase  response  of  the  sensor  was  also  experi¬ 
mentally  determined  for  an  overheat  ratio  of  0.92  and  a 
mean  shear  stress  of  30  mPa  and  is  shown  in  Fig.  14. 
The  cross-spectrum  between  the  sensor  output  voltage 
and  the  applied  acoustic  pressure  perturbation  was 
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Fig.  11  The  response  of  the  compression  driver/PWT  system  to  a 
constant  voltage  excitation  and  its  correction  (measurement  noise 
floor  included) 


Fig.  12  Gain  factor  of  the  sensor  frequency  response  function  in 
response  to  a  105  dB  (ref  20  pPa)  sine  sweep  at  an  overheat  ratio  of 
0.92  and  a  mean  shear  stress  of  50  mPa 
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Fig.  13  Linear  behavior  of  the  sensor  versus  oscillating  shear  stress 
for  a  constant  overheat  ratio  of  0.81  and  a  mean  shear  stress  of 
0.03  Pa 


taken  in  order  to  estimate  the  phase  of  the  two  signals. 
This  follows  from  the  fact  that  for  a  linear  single-input/ 
single-output  system,  the  frequency  response  function  is 
(Bendat  and  Piersol  2000) 

H(f)=W)'  (23) 

where  Gxy  is  the  cross-spectrum  and  Gxx  is  the  auto¬ 
spectrum.  Note  that  the  phase  of  the  sensor  was  not 
corrected  using  the  theoretical  relationship  between 
pressure  and  shear  stress  fluctuations.  Therefore,  at  low 
frequencies,  there  is  a  net  45°  phase  lead  as  predicted  in 
Eq.  11.  The  sensor  shows  no  resonant  peak  and  displays 
a  net  180°  phase  shift,  which  is  characteristic  of  either  a 
highly  damped  second-order  system  or  a  second-order 
system  with  two  real  poles  (i.e.,  two  time  constants).  The 
two  time  constants  can  be  explained  as  follows.  At  low 
frequencies  the  time  response  of  the  sensor  will  depend 
on  the  conduction  to  the  substrate,  but  at  high  fre¬ 
quencies  the  response  will  only  be  a  function  of  the  heat 
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Fig.  14  Phase  of  the  sensor  frequency  response  function  in 
response  to  a  105  dB  (ref  20  pPa)  sine  sweep  at  an  overheat  ratio 
of  0.92  and  a  mean  shear  stress  of  30  mPa 


transfer  associated  with  the  film  itself  (Haritonidis  1989). 
Measurements  obtained  beyond  4  kHz  display  consid¬ 
erable  scatter.  This  is  due  to  the  fact  that  the  signal  level 
is  near  the  noise  floor  of  the  uncompensated  sensor. 


5  Conclusions  and  future  work 


A  method  for  the  direct,  in  situ,  static  and  dynamic 
calibration  of  thermal  shear-stress  sensors  has  been  de¬ 
scribed.  This  technique  uses  traveling  planar  acoustic 
waves  superimposed  on  top  of  an  established  laminar 
mean  duct  flow  to  generate  a  known  oscillating  shear 
stress.  This  technique  is  novel  in  that  it  permits  inde¬ 
pendent  variation  in  the  mean  and  fluctuating  shear 
stresses.  The  ability  to  vary  these  parameters  indepen¬ 
dently  is  crucial  to  determine  the  dynamic  response  for 
various  mean  shear  stress  levels.  The  theoretical  and 
practical  aspects  of  this  technique  were  presented  as 
were  representative  dynamic  calibration  data  for  a 
micromachined  thermal  shear-stress  sensor. 

The  PWT  is  capable  of  generating  known  sinusoidal 
shear  stress  over  a  frequency  range  of  2  mPa  at  200  Hz 
to  12  mPa  at  7  kHz.  The  limiting  factors  in  the  perfor¬ 
mance  of  the  PWT  include  acoustic  reflections  especially 
at  low-frequencies  due  to  non-ideal  termination,  giving 
rise  to  a  small  percentage  of  standing  waves  and  the  flow 
transition  to  turbulence  at  low  mean  shear  stress  levels 
(~0.08  Pa).  The  presence  of  standing  waves  could,  in 
theory,  be  accounted  for  by  using  the  two-microphone 
method  (Jones  and  Stiede  1997)  to  decompose  the  field 
into  incident  and  reflected  waves.  In  addition,  the  tem¬ 
perature  fluctuation  may  contribute  to  the  dynamic  re¬ 
sponse  of  the  sensor.  Without  knowing  the  dynamic 
voltage  temperature  sensitivity  of  the  device,  a  conser¬ 
vative  estimate  of  the  voltage  fluctuation  due  to  tem¬ 
perature  can  be  estimated  via  Vj  =  %Tr,  where 


V'  **  [IbRol\ 


[(7  —  1)7- 

nr2 


(24) 


The  first  term  is  the  static  temperature  sensitivity  for  a 
fixed  bias  current  IB ,  mean  sensor  resistance  R  at  a  mean 
temperature  T,  and  a  is  the  thermal  coefficient  of  resis¬ 
tance  (Cain  et  al.  2000;  Sheplak  et  al.  2002).  The  second 
term  is  the  isentropic  temperature  fluctuation  associated 
with  an  acoustic  perturbation,  and  y  is  the  ratio  of 
specific  heats  (Morse  and  Ingard  1968).  As  discussed 
previously,  V[  =  (dF/dT)i'wall.  For  the  conditions  in 
Fig.  12,  the  ratio  V[jVj  varies  from  4  at  200  Hz  to  25  at 
7  kHz.  This  is  a  worst-case  estimate  that  will  only  in¬ 
crease  when  the  thermal  boundary  layer  is  accounted 
for.  In  any  case,  a  clear  understanding  of  the  tempera¬ 
ture  sensitivity  of  these  devices  is  required  to  accurately 
account  for  this  undesired  input. 

The  work  described  in  this  paper  is  part  of  an 
ongoing  investigation  into  the  theoretical  and  experi¬ 
mental  aspects  of  thermal  shear-stress  sensors.  The 
extension  of  thermal  sensors  towards  quantitative 
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turbulence  measurements  is  non-trivial  due  to  the 
complex  conjugated  heat-transfer  interactions  between 
the  sensor  and  the  supporting  substrate/membrane. 
In  determining  the  device  sensitivity  to  shear  stress 
fluctuations,  one  requires  an  accurate  model  of  the  fre¬ 
quency-dependent  conduction  term.  A  model  based  on 
one-dimensional  heat  conduction  into  a  semi-infinite 
medium  is  not  an  accurate  representation  of  this  term 
for  the  case  of  a  sensor  with  a  sealed  vacuum  cavity. 
Clearly,  detailed  time-resolved  flow-field/substrate  heat 
transfer  numerical  simulations  would  provide  insight 
into  the  sensor  response  to  fluctuating  shear  stress  and 
temperature. 
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